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Ultra-broadband, efficient and unidirectional surface plasmon polariton (SPP) launching is of great concern 
in plasmonic devices and circuits. To address this challenge, a novel method adopting deep-subwavelength 
slits of subwavelength period (Xspp/4 ^ ^spp/3) in a thick metal film and under backside illumination is 
proposed. A new band pattern featuring broadband and wide angular characteristics, which is due to the 
coupling of the zeroth-order SPP resonance at the superstrate-metal interface and the first-order SPP 
resonance at the metal-substrate interface, is observed for the first time in the dispersion diagram. 
Unidirectional SPP launching efficiency of ~50%, ultra-broad bandwidth of up to 780 nm, covering the 
entire optical fiber communication bands, and relatively wide angular range of 7° are achieved. This 
remarkable efficient, ultra-broadband and wide angular performance is demonstrated by carefully designed 
experiments in the near infrared regime, showing good agreement with numerical results. 

As a promising candidate to provide the new generation information network with improved bandwidth 
and speedy plasmonics has generated increasing research enthusiasms. The last decade has witnessed 
remarkable progress in the development of plasmonics and its applications in areas such as light routing 
and manipulation at the nanoscale^"'*, especially the wavelength division multiplexing (WDM) plasmonic circuit 
elements, such as the plasmonic demultiplexers^'^, broadband modulators^, and broadband antennas^. 

In order to investigate the WDM plasmonic circuits, there is a prerequisite that broadband surface plasmon 
polaritons (SPPs) propagating in a well-defined direction should be efficiently launched with appropriate 
approaches. Although various configurations have been proposed and continuous efforts have been exerted to 
establish the theory and optimize the structure, it remains challenging to realize unidirectional, efficient and 
meanwhile broadband SPP launching. The conventional prism approach making use of attenuated total 
reflection (ATR), known as the Kretschmann-Raether configuration, is very sensitive to the wavelength and 
the angle of incidence^, making it widely used in sensors. Conversely, SPP launching approaches based on one or 
two subwavelength features (such as slits, grooves or ridges) in or on a thick metal film^°"^^ may be angular 
insensitive or even relatively broadband. However, as pointed out in refs. 12 and 18, the launching efficiency 
(usually not reported) is rather small. It has been widely accepted that efficient and practical SPP launchers should 
employ periodic or aperiodic arrays of scattering objects^^. For conventional periodic gratings, the unidirectional 
SPP launching efficiency is limited although continuous efforts have been put on the geometrical optimiza- 
tion^^"^^. Recently, the depths, the widths, and the neighboring distances of aperiodic grooves have been so 
carefully optimized that —50% launching efficiency in the desired direction has been theoretically^'*'^^ and 
experimentally^^ achieved. With aperiodic metallic grooves, a record unidirectional extinction ratio of 55 has 
also been reported^^, but the coupling efficiency was not provided. More recently, the broadband and unidir- 
ectional excitation of SPPs has been achieved using chirped plasmonic gratings^ however the launching efficiency 
(not reported) is supposedly small compared with their periodic counterparts. In other words, the SPP launching 
in these periodic or aperiodic metallic gratings suffers from very limited spectral bandwidth (only a few dozens of 
nanometers) in case of high efficiency, or very low coupling efficiency in case of broad bandwidth, restricting their 
applications in WDM plasmonic circuits. A probable reason is that the periods or neighboring distances were 
usually restricted to be comparable with or larger than the SPP wavelengths. 
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In order to achieve high unidirectional efficiency and simulta- 
neously ultra-broad bandwidth of more than 415 nm so as to cover 
the entire optical fiber communication bands (O, E, S, C, L, and U 
bands), in this work, we propose a novel SPP launching method using 
periodic deep-subwavelength slits of subwavelength period (}tspp/4 
~ where Xspp is the SPP wavelength. The slits are milled in 

an optically thick metal film and obliquely illuminated from the back 
side, prohibiting the incident light from becoming significant source 
of noise and reducing the system's size^°. Interestingly, for the first 
time we find a new band pattern featuring broadband and wide 
angular characteristics, which is due to the coupling of the zeroth- 
order SPP resonance at the superstrate-metal interface and the 
first-order SPP resonance at the metal-substrate interface. The cor- 
responding unidirectional SPP launching efficiency is found to be 
—50% over an ultra-broadband range (up to 780 nm) and relatively 
large angular range (~7°). This remarkable performance makes the 
proposed SPP launching approach distinct from conventional peri- 
odic or aperiodic metallic gratings, of which periods or neighboring 
distances are usually comparable with or larger than the SPP wave- 
lengths, featuring spectrally and angularly sensitive SPP launching 
efficiencies It is experimentally demonstrated in the near-infra- 
red regime. 

Results 

The configuration in Fig. 1(a) illustrates N periodic slits of deep- 
subwavelength width w and subwavelength period p in an optic- 
ally-thick metal film of thickness h. The slits are back-side illumi- 
nated by a transverse-magnetic polarized light of oblique incidence 
angle 0. This structure is distinct from the Kretschmann-Raether 
configuration^ because of the slit array and the optically thick metal 
film. It also differs from conventional metallic gratings because p is 
much smaller than the SPP wavelength. 

We firstly calculate the band structure of the corresponding one- 
dimensional (ID) plasmonic crystal, i.e. the fully-periodic {N = oo) 
grating. To cover entire optical fiber communication bands 
(1260 nm ~ 1675 nm), we setp = 440 nm, which is about ?ispp/4 
~ }^spp/3, w = 100 nm, and h = 300 nm (optically thick enough). 
Fig. 1(b) shows the (D-k\ \ dispersion diagram, where cd is the angular 
frequency, k\ \ = ko/isubSin^ is the in-plane wavevector of the incident 
light from the substrate. The color scale represents the absorption 



intensity. As have been pointed out in refs. 27 and 28, there are 
narrow bright bands that generally follow two sets of curves: the 
folded dispersion relations of SPPs on flat metallic surfaces (black 
lines) and the Rayleigh anomaly (RA) wavelengths (red lines) at the 
Au-substrate interface (solid lines) and the superstrate-Au interface 
(dashed lines). The dispersion relations of SPPs are determined by 
the Bragg coupling condition^^. 



/coRe(nspp,sup/sub)p - koHsnhp^ine = 2mi7i, 



(1) 



where nspp,sub/sup = y SmSsub/sup/ (sm + Ssub/sup) is the effective index 
of SPP on flat metallic surface with Ssub/sup being the relative permit- 
tivity of the substrate (sub) or superstrate (sup), and mi is the dif- 
fraction order. The RA wavelengths are defined by the passing-off of 
a diffraction order (m2). 



^oSsub/supP - konsnhpsinO = 2m27i. 



(2) 



We should note that these narrow bands are observed for wave- 
lengths comparable to the periods and the SPP and RA curves almost 
merge at small energies, for which /cspp ~ /cq. These properties agrees 
well with previous studies on conventional metallic gratings^^'^^, of 
which periods are comparable with or larger than the SPP wave- 
lengths. These narrow bands explain why the SPP launching efficien- 
cies of conventional metallic gratings are spectrally and angularly 
sensitive. 

Interestingly, for the first time we observe a new bright band 
featuring large spectral and angular ranges. It locates between the 
0*^-order SPP resonance wavelength at the superstrate-Au interface, 
SPPo,sup» ^nd the first-order SPP resonance wavelength at the Au- 
substrate interface, SPPi^s^b- Its spectral and angular characteristic 
and location make it distinct from the above-mentioned narrow 
bands locates between the SPP and RA curves, as well as the slit 
Fabry-Perot (F-P) resonance curve determined by 2 ko Re(/is,eff) ^ 
+ arg(rs,top) + arg(rs^top) = 2n mi, where /is,eff is the effective index 
of the slit mode, Ts^top ^nd Ts^top ^re the reflectance coefficients of the 
slit mode at the top and bottom openings, respectively, the function 
"arg" refer to the argument of a complex number, and is an 
integer. 

To further understand the physical origin of the new band, we 
theoretically calculate the SPP launching efficiency for N periodic 
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Figure 1 | (a) Schematic of the proposed SPP launching configuration, (b) Absorption diagram of the plasmonic crystal, (c) The spectral and angular 
dependence of rj^ for N = 25 with the white contour indicating half the optimal efficiency, (d)-(f) Near-field IH^I under various wavelengths and 
incidence angles as indicated by the diamonds in (c) . The incident plane indicated by a dashed line locates at 1 00 nm below the quartz-metal interface. The 
dashed and solid arrows show incident direction and specular reflection direction, respectively. The calculations were performed with h = 300 nm, w = 
100 nm, and p = 440 nm. 
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Figure 2 | (a) Absorption diagram of the plasmonic crystal, (b) Spectral and angular dependence of rj^ for N = 25. The calculations were performed with 
h = 200 nm, w = 30 nm, and p = 290 nm. 



metallic slits, rj^ , which is defined as the ratio of the power of the 
launched forward- ("+") or backward-propagating ") SPPs to 
that of the incident light over the whole grating length, (N-l)p + w. 
Fig. 1(c) shows the spectral and angular dependence of f]^ ofN= 25. 
It reveals that f]^ reaches as large as 50% and holds more than 25% 
over extremely wide spectral range and relatively large angular 
ranges. An ultra-broad spectral full width at half maximum 
(FWHM) of up to 780 nm (from 1210 nm to 1990 nm) at ^ = 44° 
and an angular FWHM of 7° at 2 = 1550 nm are achieved. The ultra- 
broad bandwidth covers the entire optical fiber communication 
bands, making the proposed SPP launching approach very prom- 
ising in WDM plasmonic circuits. We note that rj~ is extremely 
small over the corresponding spectral and angular ranges, corres- 
ponding to large extinction ratio, which is defined as E^ = rj^ / rj~ 
and more than 20 (see Fig. SI in Supplementary Information). We 
also note that the spectral FWHM ranges of f]^ have been indicated 
by the location of the new band of the dispersion diagram in Fig. 1 (b). 

Given that the slit period of ~ >^spp/3, if and only if mi = 0, 

Eq. (1) can be satisfied. In other words, only the O^'^-order SPP res- 
onance exists for the transmitted diffraction. Each slit acts as a chan- 
nel that converts the back-side incident light into SPPs propagating 
along the superstrate-metal interface (referred to as top-SPPs). SPPs 
propagating along the metal-air interface (referred to as bottom- 
SPPs) are also generated by the slits. Meanwhile, the launched top- 
and bottom-SPPs couple with each other through the slits and also 
re-radiate to the substrate because of Rayleigh-like scattering by the 
slit array^°. These re -radiated emissions interfere with the reflected 
diffraction field by the subwavelength metallic ridges, resulting in 
energy redistribution between the reflected fleld and launched top- 
SPPs in the transmitted fleld. In the dispersion diagram, the coupling 
of top-SPPs and bottom-SPPs will result in a band between their 
corresponding dispersion curves. In other words, the new band that 
locates between SPPo,sup and SPP sub in Fig. 1(b) should originate 
from the coupling of top-SPPs at the 0^^-order resonance and bot- 
tom-SPPs at the first-order resonance. On the other hand, the inter- 
ference between SPP re-radiations and the metallic ridges' reflected 
diffraction is vividly visualized by the near fleld distributions under 
various incident wavelengths and angles, as illustrated in Figs. 1(d)- 
1(f). The incident plane locates at 100 nm below the quartz-metal 
interface so as to manifest the reflected field. It is clearly shown that 
forward-propagating SPPs are launched efficiently while backward- 
propagating SPPs are very weak over the entire optical fiber com- 
munication bands and at various incident angles. 

Since the new band locates between SPPo,sup and SPPi^sub^ it is 
possible to tune its location and accordingly the ultra-broad band- 
width of f]^ to the desired spectral range by varying SPPo,sup and 
SPP 1, sub J which are determined by structural parameters, such as p, 
^sup and Wsub- By decreasing the geometrical parameters (p, w, and h), 



the location of the new band is shifted to higher frequency, and 
accordingly the bandwidth of rj^ is blue-shifted. Fig. 2 shows the 
results for h = 200 nm, w = 30 nm, and p = 290 nm. Comparing 
Figs. 1(c) and 2(a), we observe that the spectral bandwidth has been 
blue-shifted from 1210 nm ~ 1990 nm to 890 nm ~ 1270 nm at 
0 = 44° by reducing the slit period from p = 440 nm into p = 
290 nm. It is worth to note that the slit period holds to be about 
/lspp/4 ~ /lspp/3. The blue-shift by decreasing geometrical size is in 
accordance with the scaling law^\ Moreover, the angular FWHM is 
about 7° at >^ = 1000 nm. 

Note that the idea of the proposed SPP launcher originates from 
the dispersion properties of an inflnite periodic array of slits, while 
the use of the SPP launcher should adopt a flnite number N of slits so 
as to excite unidirectional SPPs on flat metallic surfaces. For conven- 
tional SPP launchers composed of periodic features with period 
comparable with >^spp» the SPP launching efflciencies depend on 
A^l This property also holds for the proposed SPP launcher with 
some differences. As shown in Fig. 3(a), f]^ flrst increases to about 
51% and then decreases slowly as AT increases. The decrease of f]^ is 
because the saturation of the launched SPPs due to scattering and 
propagation losses. Compared with conventional SPP launchers, 
which exhibit low f]^ and a rapid saturation of f]^ as N increases^^ 
the proposed SPP launching configuration is of much higher effi- 
ciency and much slower saturation. Fig. 3(a) also shows that the 
extinction ratio is larger than 100 if AT > 6, indicating that the pro- 
posed configuration launches highly unidirectional SPPs. 
Comparing Figs. 1(c), 3(b) and 3(c), we notice that although the 
angular FWHM of the proposed SPP launcher decreases and the 
optimal efficiency varies as N increases, the ultra-broadband per- 
formance holds for various N. In other words, the efficient, unidir- 
ectional, ultra-broadband and relatively wide angular performance 
holds for the proposed SPP launchers of a very wide range of N. For 
the design of a SPP launcher of given AT, the desired spectral ranges 
and the corresponding/), one may optimize w and h for the optimal 
efficiency using efficient theoretical models introduced in ref. 32. 

To experimentally demonstrate the remarkable performance, we 
carefully designed samples composed of a launcher and a decoupler, 
as schematically shown in Fig. 4(a). Under back- side illumination of 
a TM-polarized light beam, forward-propagating top-SPPs are 
excited with a launching efficiency r]^ , and then scattered into far 
field by the decoupler with a decoupling efficiency rj^. The decoup- 
ling efficiency r]^ is defined as the ratio of the re- radiated power Pq 
from the decoupler within a cone with a ±30° extraction angle, 
which is determined by the numerical aperture (N A) of the collecting 
microscope objective (NA = 0.5 here), to the power of top-SPPs at 
the slits. In Fig. 4(a), Pin is the incident power onto the launcher. 

Taking into account of the tunable wavelength range of the laser 
source we used in the experiment (800 nm to 920 nm) and the 
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fabrication yield, we designed the launcher of a sample, Launcher I, 
to be of the following parameters: w = 60 nm,p = 240 nm, and N = 
15. The launcher of the other sample. Launcher II, has the following 
design parameters: w = 120 nm, p = 480 nm, and N = 8. Both 
launchers have the same slit length (L = 54 |im). Furthermore, the 
decouplers share the same parameters as Launcher II, except for the 
slit length L/2, because of two considerations: the noise power could 
be eliminated and the SPP coupling efficiency could be accurately 
inferred, as will be elaborated later. The scanning electron micro- 
scopy (SEM) images are shown in Figs. 4(b) to 4(d) (and Figs. S2(a)- 
(b) in Supplementary Information). From the calibrated SEM 
images, the measured slit parameters are determined with the aver- 
age values and the standard derivations, which are obtained from 
statistics of slit widths and distances at different positions: w = 55 ± 
15 nm andp = 248 ± 15 nm for Launcher I, w = 105 ± 6 nm andp 
= 495 ± 7 nm for Launcher II, and w = 120 ± 7 nm andp = 499 ± 
7 nm for the decouplers. 

With the optical setup schematically illustrated in Fig. 4(e), we 
obtained the images of the decoupler and Launcher I with a charge 
coupled device (CCD), as shown by Figs. 4(f) and 4(g), respectively. 
We note that the far-field reflection by Launcher I is significantly 
suppressed compared with the specular reflection by the surround- 
ing flat metal surface, as predicted by the near-field distributions 
shown in Figs. 1(d) -1(f), although different parameters were used. 
The CCD images were further processed (see Methods) to accurately 
obtain the total launching- decoupling efficiency expressed as 
rj^ =^^yj^ exp ( — ///spp)j where /spp is the SPP propagation length. 
We measured /spp adopting a method improved from the one intro- 
duced in ref 33. Fig. 5 shows the structure design, a fabricated sam- 
ple, the optical characterization and the image post-processing (see 
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Figure 3 | (a) Dependence of rj^ and Er on N. (b)(c) Spectral and angular 
dependence of rj^ for AT = 5 (b) and iV = 50 (c). The calculations were 
performed with /I = 1550 nm, 9 = 44°, and other parameters (h, w, p) are 
the same as Fig. 1(c). 



Methods), and the measured results comparing with the theoretical 
values calculated using /spp = 1/(2 /cq lm{n^^^)). The shorter- than- 
expected /spp in Fig. 5(d) is probably due to inelastic SPP scattering by 
the surface roughness inherently introduced during evaporated and 
increased by FIB patterning^^. To infer f\'^ accurately from the mea- 
sured =^c^d ^^p( ~ Vfepp)' Launcher II and the decouplers are 
carefully designed so that an almost constant relationship exists 
between rj^x and t]^ (11,6 = 45°) y where II (or I) represents Launcher 
II (or Launcher I), over the laser wavelength range in our experiment, 
as shown in Fig. 6(a): 1 .19 f]^ (11,^ = 45°). The linear relationship 
holds for the fabricated samples, but it is slightly modified into 
rj^^l.lVrj^ (11,0 = 45°). As a result, we assume that this approxi- 
mation also holds for the measured efficiencies and then infer rj^ 
and f]d from the measured f]^ as: ^^^^—^1.17?/+ (11,^ = 45°) 
exp(-///spp)and (I/II)-f^+ (I/II)/^/^ exp(///spp). 

Figures 6(b) and 6(c) compare the experimental and the calculated 
T]^ of the two launchers with respect to IforO = 45° and to ^ for >^ = 
830 nm, respectively. Although the measured efficiencies are smaller 
than the calculated results, especially for Launcher I, they agree well 
in general trends. In Fig. 6(b), a slight oscillation occurs for the 
experimental data of Launcher I because of the interference between 
top-SPPs and bottom-SPPs at the decoupler, which should not 
appear in practical applications. Limited by the wavelength range 
of the laser source we used, we have only demonstrated the front 
part of the broadband and meanwhile efficient performance (see 
Supplementary Information) in the near infrared regime. In 
Fig. 6(c), the angular FWHM of more than 8° has been demonstrated 
experimentally. 

Discussion 

The difference between measured and calculated efficiencies in 
Figs. 6(b) and 6(c) could originate from the fabrication imperfec- 
tions. As we have pointed previously, the measured slit parameters 
determined from the calibrated SEM images indicate that the fab- 
rication deviations of Launcher I are relatively large. This is due to the 
large size of the whole sample and the deep -sub wavelength slit width. 
Aside from the size deviations, some metallic ridges of Launcher I 
separate from the substrate because of extremely small widths (see 
Fig. 5(d)). However, even with these fabrication imperfections, all the 
measured rj^ (I) are larger than the measured slit duty factor (22.2%) 
across the laser wavelength range in the experiment for ^ = 45°, and 
the maximum value is up to 31.9%. 

Up to now, the largest efficiency of —52% has been reported from 
a carefully optimized aperiodic SPP launcher^^, which is illuminated 
from the top and sensitive to both 1 and 0. As pointed out in ref 10, 
the incident light from the top is a significant source of noise, unless 
directed away from a region of interest, which then decreases the 
signal and increases the system's size. This problem is eliminated 
using back-side illumination in optically thick metal films. 
Moreover, the neighboring distances of the aperiodic SPP launcher 
in ref 18 is comparable to the SPP wavelength, whereas the period in 
our proposed method is of /lspp/4 ~ >^spp/3. As a result, the proposed 
SPP launching approach adopting back- side illumination and sub- 
wavelength period will also gain advantages in integration besides the 
attractive ultra-broadband and meanwhile efficient performance. 
Additionally, compared with aperiodic gratings^' ^^'^^'^^ of which the 
widths, the depths and the neighboring distances are different, the 
proposed periodic sUts are much easier to fabricate. 

In summary, we have proposed and demonstrated a novel method 
to realize unidirectional, efficient and meanwhile ultra-broadband 
SPP launching adopting deep -sub wavelength metallic slits with per- 
iod of >^spp/4 ~ >lspp/3. We expect that the proposed SPP launching 
method will be favorable in nanophotonic applications, especially 
WDM plasmonic circuits. 
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Figure 4 | (a) Schematic of the sample layout, (b)-(d) SEM images and close-ups of a sample composed of the decoupler and Launcher L (e) Schematic of 
the optical characterization setup, (f)-(g) CCD images of the decoupler and Launcher I, respectively. 



Methods 

Simulation. The band structure of the ID plasmonic crystal, the SPP launching 
efficiency and decoupling efficiency are calculated using the rigorous coupled 
wave analysis/aperiodic Fourier modal method (RCWA/a-FMM)^*■^^ The near fields 
are calculated using Lumerical FDTD Solutions. The wavelength-dependent complex 
permittivities of gold are interpolated from experimental data^^, the superstrate and 
the material filling the slits are assumed to be air (/^sup — l-0)> «^nd the refractive index 
of the substrate n^^\, = 1.46 is used throughout this work. 

Fabrication. A 200 nm-thick gold film is sputtered onto a 2 cm X 2 cm quartz 
substrate. A focused ion beam (FIB) milling system is used to mill samples composed 



of a launcher and a decoupler with a separation distance /, and samples composed of a 
central slit and two surrounding grooves, all of which are 120 nm wide and separated 
with various pairs of asymmetric slit-groove distances (9, 12, 15, 21, 24, and 27 |am). 

Characterization. The optical measurement is conducted using home-built 
experimental setups to obtain the CCD images. The light source is a continuous wave 
Ti:sapphire laser with a tunable wavelength range of 800 nm to 920 nm. The laser 
beam expanded by a telescope system and attenuated by a calibrated vari-ND filter 
illuminates a rectangle iris, which is imaged onto the launcher or the central slit from 
the substrate side using an imaging system composed of a lens Li and a microscope 
objective Oi (SOX, NA = 0.5). To measure the SPP coupling efficiency, the sample is 




800 820 840 860 880 900 920 
X (nm) 

Figure 5 | (a) SEM image of a slit-groove doublet sample with di= \2 |am and d2 = 24 \im. (b) (c) CCD images of the sample without and with the spatial 
filtering system, respectively, (d) Theoretical and measured /spp versus X with error bars denoting the standard error of the mean (SE). 
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Figure 6 | (a) Numerical spectral behaviors of r]^ (11,6 = 45°) and r]^ 
calculated for the design parameters (w = 120 nm and p = 480 nm for 
both Launcher II and the decoupler) and the average measured parameters 
( w = 105 nm and p = 495 nm for Launcher II, and w = 120 nm and p = 
499 nm for the decoupler. (b)(c) Experimental and a-FMM calculated rj^ 
versus /I at 0 = 45° (b), and versus 0 at /I = 830 nm (c). The a-FMM 
calculations were performed with the average of measured slit parameters. 



attached to a prism with index matching oil and then mounted on an XYZ linear 
translation plus Z axis rotation stage, so that the angle of incidence could be varied 
continuously. The decoupler, the launcher, and the iris's image reflected by the flat 
gold film are separately imaged onto a CCD using another imaging system composed 
of an objective O2 (50 X , NA = 0.5) and a lens L2. Another calibrated vari-ND filter is 
placed in front of the CCD to keep the CCD away from saturation and to improve the 
signal-to-noise ratio (SNR) of obtained CCD images. To further improve the SNR, a 
spatial filter is used to block the reflection of the iris's image when the decoupler is 
imaged. This is because the re-radiation of the decoupler is very weak compared with 
the reflection of the iris's image without the spatial filter (see Fig. S2(c) in 
Supplementary Information). To measure the SPP propagation length, the sample is 
freestanding mounted. Normally incident light passed through the central slit 
launches bidirectional top-SPPs of the same intensities, which are scattered by the 
surrounding grooves (see inset in Fig. 5(a)). The re-radiations of both grooves are 
collected by another imaging system and imaged onto the CCD, whereas light that 
directly transmits through the slit is blocked by a spatial filtering system to improve 
the SNR. 

CCD image processing. For the measurement of SPP launching efficiencies, the CCD 
images are processed as follows. The raw power of the re-radiation from the decoupler 
Pr is calculated by integrating the pixel intensities within the blue dashed box in 
Fig. 4(f), which is occupied by the decoupler and obtained via image edge recognition. 
The power of noise is determined by adding the pixel intensities of the two yellow 
dashed boxes in Fig. 4(f). The yellow boxes located at the two sides of the decoupler 
are of the equal width but of half the length of the blue box. The noise-free radiation 
power is then obtained as Pq = Pr - Pn- The incident power Pin is obtained by 
measuring the reflection of the incident beam from flat gold film, i.e. integrating the 
pixel intensities in the same area as the white dashed box occupied by the launcher in 
Fig. 4(g). The white box occupied by the launcher is also obtained via image edge 
recognition. Note that the calculation of Pq and Pin have taken into account the fact 
that CCD images of the decoupler and Launcher I are obtained using different 
transmission efficiencies of the calibrated vari-ND filter in front of the CCD, Tq and 
Tin, respectively. Taking into account that the decoupler is of half the length of the 
launcher, the total launching- decoupling efficiency incorporating the SPP 
propagation loss is written as rj^ =^c^d ^^P (~ V^spp) = 2Po/Pin- The image 
processing for measuring the propagation length is similar. In Fig. 5(c), Pri and Pr2 are 
raw powers of the re-radiations from the grooves indicated by blue dashed boxes, 
which are extracted via image edge recognition. The noises powers P^i and Pn2 are 
determined by adding the pixel intensities of the yellow- dashed boxes. The yellow 
boxes located at the two sides of the blue ones are of equal widths but are half the 
length of the blue boxes. As a result, Pgi = Pri - Pni and Pg2 = Pr2 - Pni are the 
respective noise-free re-radiation powers from the left and right grooves. At every 



wavelength, we fit the decay of the re-radiation power with an exponential to extract 
the top-SPP propagation length, i.e., /spp = {d2 - (ii)/(lnPgi - lnPg2). 
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